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For systems of functional differential equations that can take into 
account finite or infinite delays, a complete characterization is given 
of the invariance of positivity in the sense that all solution components 
are positive whenever the initial condition function is positive. A 
related result concerning a comparison of the solutions of pairs of 
initial value problems is also given. One application of the results 
described concerns a model for synchronizing geographically sepa- 
rated oscillators and another is in the area of economics. 

I. INTRODUCTION 

Consider a system of functional differential equations of the form 

x = f(t,x t ), t>to, *,„ = <*>, (1) 

in which x is a real n-vector valued function of t, x denotes dx/dt, 4> is 
an initial condition function, and x, denotes the function defined by 
Xt(s) = x{t + s) for s < 0.* (When f(t, x,) depends only on t and x,(0), 
(1) reduces to a system of ordinary differential equations.) 

The main purpose of this paper is to give a solution to the problem 
of determining conditions under which (under certain typically very 
reasonable conditions on f), x{t) of (1) has components that are all 
positive for t > to whenever <J> is positive in, for example, the sense that 
<£(s) has positive components for s < 0. The problem arises in connec- 
tion with the mathematical modeling and analysis of economic proc- 
esses, and it comes up in several other areas as well. (An example 
concerning the synchronization of geographically separated oscillators 
is described in Section 2.6.) In some instances, the invariance of 
positivity in the sense described above is crucial, in that the lack of 

* For background material concerning equations of the type (1), see, for example, 
Refs. 1 and 2. 
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positivity of a component of x(t) for some t and positive <£ means that 
the associated model is inappropriate. 

Our main result, Theorem 1 of Section II, is concerned with the case 
in which / is continuous and locally Lipschitz. It provides an explicit 
and useful condition under which positivity is invariant, and it also 
asserts that positivity is invariant if and only if (1) preserves non- 
negativity in the sense that x(t) has nonnegative components for t > 
to whenever <£ is nonnegative. 

The nonnegativity-preservation problem has been considered in 
Refs. 3, 4, and 5, and the relationship between Theorem 1 and the 
earlier material is indicated in Section 2.3. 

A corollary of Theorem 1 is as follows. Suppose that (1) is a system 
of ordinary differential equations, and that /is continuous and satisfies 
a global Lipschitz condition (in the usual sense). Let g(t, x) denote 
f(t, x t ). Then positivity is invariant for (1), by which we mean invariant 
for each starting point to, if and only if for each i, gi(to, v) > for each 
to and each real n-vector v such that Vi = and vj > for/ # i. Notice 
that for the special case in which f(t, x,) = Ax, where A is an n x n 
matrix of real constants, our condition is equivalent to the requirement 
that the off-diagonal elements of A are nonnegative. The corresponding 
proposition concerning nonnegativity preservation for this case is well 
known. 6 Of some interest is the fact that the corollary described above 
becomes false if the Lipschitz hypothesis is replaced with the assump- 
tion that (1) has exactly one solution for each initial condition (see 
Section 2.3). 

A result related to Theorem 1 that provides a necessary and suffi- 
cient condition for the invariance of positivity, or of nonnegativity, of 
the difference of the solutions of a pair of equations of the type (1) is 
given in Section 2.4. Specific applications of that result, as well as of 
Theorem 1, are described in Section 2.6. 

II. CHARACTERIZATION OF THE INVARIANCE OF POSITIVITY 
2.1 Preliminaries 

We use the following notation and definitions. With n an arbitrary 
positive integer, R" denotes the set of real n-vectors v = (v\, Vz, •••, 
v n ), Rl = {v e R": Vi > for each i), and \v\ = max, 1 1>, | for v e R". For 
u and v in R", the inequality u > v (u > u) means that m, > u, (w, > d,) 
for each i. The zero n-vector of R" is denoted by 8. 

We denote by C the Banach space of bounded continuous functions 
from (—oo, 0] to R", with norm given by 

| w | = sup | w(t) | 

I f (-», 0| 

for all w e C. 
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The symbol T denotes any real interval of the form [a, oo), (a, oo), or 
(—oo, oo), and to is an element of T. For each t e T and each bounded 
continuous function w from (-00, t] to R", w, denotes the element of C 
defined by w t (s) = w(t + s) for s < 0. 

Throughout Section II, /"in (1) denotes a mapping of T X C into R". 
We say that f is continuous in t if f(t, w,) is a continuous function of t 
for t > to whenever t e T and w is a bounded continuous mapping of 
(—00, 00) into R", and we say that / is locally Lipschitz if for each 
t e T, each y e [to, °°), and each compact set B in R", there is a constant 
p(to, y, B) such that \f(t, u) - f(t, v)\ < p(to, y, B)\u - v\ for each 
t e [to, y] and each u and v in C such that the range of u, and also of v, 
is contained in fi. 

A solution of (1) through a given (to, <i>) e T X C means a continuous 
iT-valued function x that is defined on (-00, 00), is differentiate on (to, 
00), and is such that (1) is satisfied (with the understanding that at 
t = to, x denotes the right-hand derivative).* As in the case of ordinary 
differential equations, if f is continuous in t and satisfies a uniform 
Lipschitz condition in the sense that / satisfies a local Lipschitz 
condition with p(t , y, B) independent of B, for each (to, <j>) e T X C 
there is a unique solution of (1) through (to, <J>) (see p. 409 of Ref. 2). 

In the next section, we refer to the following hypothesis. 
H.l: There is a solution x of (1) through each (to, <j>) e T X C, and f is 
locally Lipschitz as well as continuous in t. (In particular, each 
solution of (I) is unique.) 

2.2 Our principal result 

Under the assumption that H.l holds, consider the following prop- 
erties and condition. 

Property 1 (Invariance of Positivity, Version 1): For each (to, <J>) e T 

X C such that <f>(0) > 6 and 4>(s) e Rl for s < 0, we have x(t) > for 

t>U. 

Property 2 (Invariance of Positivity, Version 2): For each (to, <i>) e T 

X C such that <f>(s) > for s < 0, we have x(t) > 6 for t> to. 

Property 3 (Invariance of Nonnegativity): We have x(t) >0fort>t Q 
whenever (to, <f>) e T X C with cj>(s) e Rl for s < 0. 

Condition 1: For each i, f(t , </>) > whenever (to, <t>) e T X C with $(s) 
ER" + fors<0 and 0,(0) = 0. 

Theorem 1: Let H.l be satisfied. Then the following four statements 
are equivalent: Property 1 holds, Property 2 holds, Property 3 holds, 
and Condition 1 is met. 



* It will become clear that our development can be extended at once to cover the case 
in which a solution need be defined on only an interval of the form (-co, /?) with ft > to. 
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Proof. We first show that Property 3 and Condition 1 are equivalent. 

Suppose that Condition 1 is satisfied, that (to, </>) e T X C is given, 
with <Ms) e Rl for s < 0, and that x(t) > for t > t Q is violated. Then 
there is a t' > fo and an index <f such that */ (£') < 0. Since /" is 
continuous in t and locally Lipschitz, the Bellman-Grownwall Lemma 7 
and Theorem 3 of Ref. 2 can be used to show that, given any r > t', 
there are an e > and an i?"-valued function w defined on (— oo, t], and 
differentiate on (to, t), such that ibi = f(t, w,) + e, t e [to, t), i = 1, 2, 
• • •, n, with w ti = <f> and w, (V) < 0* Let J = {i: Wi(t) < for some t e 
(to, t)}, and for each i e I, let U = mf{t e (to, t): wAt) < 0}. Choose k so 
that tk = min {£,: i e I}. We have W/,(t/ t ) = 0, ibi,(h) ^ 0, and 
Wf t (s) e Rl for s < 0. Thus, f tl (tk, w tl ) + e < 0, which contradicts 
Condition 1. Therefore, we have Property 3 when Condition 1 is met. 
On the other hand, if Condition 1 is not satisfied, there is an index £ 
and a (to, <f>) e T X C with <}>(s) e Rl for s < and <v (0) = such that 
fAto, (j>) < 0. Since there is a solution x through (to, $), and it clearly 
satisfies x, (t) < for (t — to) positive and sufficiently small, we see 
that Property 3 implies that Condition 1 is met. This proves the 
equivalence of Property 3 and Condition 1. 

It is clear that Property 1 implies Property 2. To see that Condition 
1 is satisfied when Property 2 holds, suppose once more that Condition 
1 is not met. Then, as in the paragraph above, there is a (to, <f>) e T X 
C with <f>(s) e R% for s < such that x satisfies x, (f) < for some <?and 
t' > to. By H.l and Lemma 2 of Ref. 2 (which is a result concerning the 
continuous dependence of the solution on the initial data), there is a 
(*o, $) e T X C such that <j>(s) > 6 for s < and such that the 
corresponding solution x meets X/ (t') < 0. Therefore, Condition 1 is 
satisfied when Property 2 holds. To complete the proof of the theorem, 
we now show that Property 3 implies Property 1. 

Assume that Property 3 holds and that Property 1 does not hold. 
Then there isa (to, $) e T X C and a corresponding solution x such 
that 0(0) > 0, <t>(s) e R" + for s < 0, x(t) eRl for t > to, and x, (f) = for 
some /and t ' > to. We assume without loss of generality that t ' — \x\i[t 
> t : Xi(t) = for some i). Thus, x(t) > d for t e [to, t'). 

Let 

B = [veRl: v; < sup\x(t)\, i = 1,2, •••,nl. 

From x,(t') = and the observation that 
"' xAs) 



I 



r x,(s) 



ds = ln[x, (t)] - ln[*, (t")l t e [t", t) 



* At t = to, the unique solution with which Theorem 3 of Ref. 2 is concerned has a 
right-hand derivative equal to the value of the functional that corresponds here to /. 
(See the proof of Theorem 2 of Ref. 2.) 
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for any t" e (to, t'), we see that x f (t) [x, (t)]~ l is not bounded from 
below on (to, t'). Let a e (to, t') be chosen so that x, (a) + p(t , t',B) 
X/-(a) < 0, in which p is the Lipschitz constant in Section 2.1. 

Define w(a) e C by [10(a)], (s) = max[0, (x„), (s) - (x a ), (0)] for s < 0, 
and [w(a)]i(s) = (x„)i(s) for s < and i ^ £. It can be verified that 
w(a)(s) e B for s < 0, and that | x - w(a) | < x, (a). 

Therefore, 

f f \a, w(a)] = x, (a) - f, (a, x„) + ft [a, w(a)] 

< x, (a) + p(t , t', B) \x„ - w(a) | 

<x, (a) + p(t ,t',B) x,(o). 

We have f,[o, w(a)] < 0, (a, w(a)) e T X C, w(a)(s) eRI for s < 0, and 
[w(o)]s(0) = 0, which, in view of the equivalence of Condition 1 and 
Property 3, contradicts our assumption that Property 3 holds. Thus 
Property 3 implies Property 1, which completes the proof of Theo- 
rem 1. 

2.3 Notes 

The Condition 1-implies-Property 3 assertion of the theorem be- 
comes false if the hypothesis that /"is locally Lipschitz is dropped and 
Property 3 is modified in the natural way so that it concerns all 
solutions that correspond to the indicated type of initial condition. 4 

The following example shows that the theorem becomes false if the 
Lipschitz hypothesis is replaced with the assumption that (1) has at 
most one solution for each (to, 4>) e T X C. Let n = 1, and let /be 
defined for all t by f(t, x,) = -(x(t)) l/2 for *(*) > 0, and f(t, x t ) = for 
x(t) < 0. Observe that /"is continuous, and that a solution (in the usual 
sense) of x = f(t, x,) for t > 0, x(0) = x°, is given by x(t) = x° for t > 
if x° < 0, and x(t) = ((x°) 1/2 - Vd) 2 for t e [0, 2(x°) 1/2 ] with x(t)=0 
for t > 2(Jt°) 1/2 if *° > 0- It can De verified that there are no other 
solutions, even though /is not locally Lipschitz. While here x(t) is 
nonnegative for t > whenever x° > 0, it is obviously not true that 
x(t) is positive for all t > whenever x° is positive. 

Essentially, the fact that Condition 1 and Property 3 are equivalent 
for ordinary differential equations is proved in Ref. 4, and in Ref. 5 
(and in the setting provided by the results in Ref. 2) that result is 
extended to cover the more general case. At the time Ref. 4 was 
written, this writer was unaware of Ref. 3, which contains a theorem 
(proved in a very different way) from which the result in Ref. 4 can be 
obtained. Our proof of the equivalence of Condition 1 and Property 3 
is basically the same as the proof in Ref. 4 for the ordinary differential 
equations case. We did not omit the proof mainly because a modifi- 
cation of it is referred to in the next section. Also, for the case in which 
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(1) takes into account only finite delays (i.e., for equations of "retarded" 
type), a direct variation of the proof, using the continuous dependence 
result in, for example, Ref. 1, p. 41, shows that Condition 1 and 
Property 3 are equivalent without the Lipschitz hypothesis, provided 
that (1) has exactly one solution for each (to, <f>) e T X C. 

Our proof of Theorem 1 shows also that, when H.l is met, Condition 
1 is necessary and sufficient that positivity is invariant in the sense 
that Property 3 holds and we have x t (t) > for t > to whenever (to, 
<f>) e T X C is such that <t>(s) e Rl for s < and the index / is such that 
<M0) > 0. 

2.4 The comparison theorem 

The proof of Theorem 1 can be modified to establish a corresponding 
theorem concerning a comparison of the solutions of two initial value 
problems. To describe that result, let g be a function from T X C into 
R", and consider together with (1) the equation 

y = g(t,y,), t>to, #„-*, (2) 

as well as the following hypothesis, properties, and condition. 

H.2: For each (to, <f>) e T X C, (1) has a solution x, and similarly, for 
each (t Q , \p) e T X C, (2) has a solution y. The mappings f and g are 
continuous in t, and at least one of the mappings f or g is locally 
Lipschitz. 

Property 4: For each (t , <f>, \p) e T X C X C such that <f>(0) > \p(0) and 
(f>(s) > \p(s) for s < 0, we have x(t) > y(t) for t > to (i.e., we have x(t) 
> y(t) for t > to for any solution x of (1) through (to, <J>) and any 
solution y of (2) through (to, $)). 

Property 5: For each (t , <f>, ^) e Tx C X C such that <j>(s) > i//(s) for s 
< 0, we have x(t) >y(t) for t > t . 

Property 6: We have x(t) >y(t) for t> to whenever (t , $, »//) e T X C 
X C such that <f>(s) > >//(s) for s < 0. 

Condition 2: For each i, f(to, <t>) > gi(to, *P) whenever (to, <j>,\p)eTxC 
X C with 4>,(0) = &(0) and <|>(s) > i//(s) for s < 0. 
Our result is the following. 

Theorem 2: If H.2 is met, then Property 4, Property 5, Property 6, and 
Condition 2 are equivalent to one another. 

Proof. The proof is similar to the one given of Theorem 1. In fact, 
straightforward modifications show that Property 6 and Condition 2 
are equivalent, and that Property 5 implies Condition 2. Since it is 
clear that Property 4 implies Property 5, it therefore suffices to use 
the equivalence of Condition 2 and Property 6 to prove that Property 
6 implies Property 4. We do that as follows. 
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Assume that Property 6 holds, but that Property 4 fails to hold. 
Thus there are (to, <f>, \p) e T X C X C, a number t' > t , and an index 
£ such that 4>(0) > ^(0), <j>{s) > +(s) for s < 0, x, (f) = y, (f), and x(t) 
>y(t) forte [to, t').Let 

B= {t; e i? n : M < sup[max ( |x(f) | , \y(t) | )]}. 

Assume for the moment that /"is locally Lipschitz. Using the fact that 
\if(t) - yAt)][xAt) - y< (t)]~ l is not bounded from below on (to, t') 
(see the proof of Theorem 1), choose a e (t , t') so that x,(o) — y,(a) 
+ p(ta, t',B)[x,(o)-y,(o)]<0. 

Define u(o) e C by the conditions [u(o)], (s) = max[(y„), (s), (x„), (s) 
- MAO) + (y„), (0)], s < 0, and [u(a)],(s) = (x a ),(s), s < 0, i * t. We 
have [u(a)]A0) = (v„),(0), and y a (s) < u(a)(s) < x a (s) for s < 0. In 
particular, u(a)(s) e B for s < 0. Also, if («„)/ (s) - (jc„), (0) + (v„), (0) 
> ( v ), (s), then (*„), (s) - [u(a)], (s) = (x„), (0) - (y„), (0). On the other 
hand, if (x„)As) - (x,)A0) + (y„M0) < (y )As), then (x )As) - (y a )As) 
< (x o )A0) - (y a )A0). Consequently, | x„ - u(a) \ < (x a )A0) - {y,)A0). 
Therefore, 

f,[o, u{a)] - g ( {a, y a ) = x t (a) - y, (o) + f, [a, u(a)] - f, (a, x ) 

< x, (a) - y ( (a) + p(t Q , t', B)[x, (a) - y, (a)], 

which shows that f f [a, u(o)] - g, (a, y„) < 0. This contradicts Condition 
2 and hence Property 6. A similar contradiction can be obtained when 
g rather than f is locally Lipschitz. (In the analogous argument, the 
function v(a) e C that plays the role of u(a) is defined by [v(o)], (s) = 
min [(x,)As), (y a ), (s) - 00/ (0) + (x a ), (0)] for s < 0, and [u(o) ],-(*) = 
(y a )i(s) for s < and i ¥> I) This shows that Property 6 implies Property 
4, and it completes the proof. 

2.5 Comments 

Since Property 6 does not hold when f = g and (1) has more than 
one solution through some (to, </>) e T X C, we see that the Condition 
2-implies-Property 6 part* of Theorem 2 becomes false if the hypoth- 
esis that at least one of the functions f or g is locally Lipschitz is 
omitted. The example given in Section 2.3 shows that the theorem 
becomes false even if the hypothesis is replaced with the assumption 
that (1) and (2) have at most one solution for each (to, <t>) e T X C and 
each (to, yp) e T X C, respectively. On the other hand, the equivalence 
of Condition 2 and Property 6 holds for equations of retarded type 
when the Lipschitz hypothesis is replaced with the assumption of 
uniqueness of solutions for at least one of the equations (1) and (2).f 



* For ordinary differential equations with f = g, this part of Theorem 2 is along the 
lines of a well-known result (Ref. 8). 

t See the corresponding comment in Section 2.3. 
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The proof of Theorem 2 described in the preceding section can be 
used to verify that C in Theorem 2 can be replaced with the set of 
continuous bounded functions from (—00, 0] to D, where D is any open 
rectangular interval of R", provided that by a solution of (1) or (2) is 
meant a solution whose values are contained in D for t > to. (In this 
connection, note that u(a) of the proof of Theorem 2 satisfies u(a)(s) 
e D for s < whenever D C R" is a rectangular interval and both x„(s) 
and y a (s) are contained in D for s < 0.) Of course, the case in which D 
= {v e R n : v > 6} is of particular interest. 

2.6 Applications 

There are many applications of Theorems 1 and 2. As a simple 
example for the purpose of illustration, consider the delay-differential 
equations 

n 

xM) = £ hij{t)[xj(t - Ty) - Xi(t)], t > 

1-1,2, ...,», (3) 

which arise 9 in the study of models of systems for synchronizing 
geographically separated oscillators. In (3), each Ty is a nonnegative 
constant, each hy is nonnegative, continuous, and bounded on [0, 00 ), 
Xj(t) denotes the frequency of the ith oscillator, and the h,j can depend 
on x as well as on certain fixed nonlinear functions. Under certain 
reasonable hypotheses concerning the hy (see Ref. 9), given a contin- 
uous x{t) for t e t, where t = [— max,*,- t,>, 0], there is a constant p such 
that for each i, Xi(t) —> p as t—> ™. Assume here that there is such a p 
for each initial-condition function. 

Theorem 1 shows that each Xi(t) in (3) is positive for t > whenever 
x(t) > 6 for t e t.* (The nature of the dependence of the ha on x is not 
of consequence at this point. If it were not true that x(t) > 9 for t > 
whenever x(t) > 6 for t e t, we would have a contradiction to the 
theorem.) Assuming now that the hij are independent of x, it follows 
from Theorem 2 that, for example, p is either increased or unchanged 
when x(t) for t e t is replaced with any continuous x(t) for which x(t) 
> x(t) for t E T. 

Two related observations concern the equations 

n 

Xi(t) = -b i0 [xi(t)] + X bij[xj(t - t v )] -I- uM), t > 
/-1 



Mi 



i=l, 2, ...,/i (4) 



* This proposition is a special case of Lemma 1 of Ref. 9 whose proof is very different. 
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of a model of a compartmental system with delays, 10 in which the 6,- 
and the by for i ^ j are locally Lipschitz monotone-nondecreasing 
functions such that 6, (0) = &„(0) = 0, each u, is continuous and 
satisfies uAt) > for t > 0, and, as in the preceding example, the Ty are 
nonnegative constants. In (4), Xi(t) denotes the amount of material in 
the ith compartment. 

From Theorem 2, we see that if x" is a solution of (4) corresponding 
to (mi, u 2 , ■-, u n ) = u" and x(t) = x a {t) for t e t, where again t = 
[-max,>, T,y, 0], and similarly with regard to x h and u b , then we have 
x a (t) > (>) x h (t) for t > when u"(t) > u b {t) for t > and x"(t) > (>) 
x b (t) for t e t. The "> part" of this proposition was given in Ref. 10. 
From Theorem 1, it is clear that we have x(t) > for t > whenever 
x(t) > for t e t, which does not seem to have been proved earlier, 
even for the case in which t,> = for all i ¥=j.* 

Consider now the case in which /"in (1) is given by f,(t, x,) = h,(x,) for 
each i, where each hi is a functional on C with the property that h,(u) 
> hi(v) for all u and v in C such that u,(s) = vfa) and u(s) > u(s) for s 
< 0. Functions /" of this form are generalizations of time-invariant 
quasimonotone (or Wazewski-type) functions" which are of interest in 
several areas, including economics. In economics applications, the Xi(t) 
often denote prices (see, for example, Refs. 12 and 13). Observe that 
(4) is of the form considered here when the u, are independent of t. 

Assuming that H.l is met, Theorem 1 provides the following simple 
necessary and sufficient condition for the invariance of positivity in 
the sense of Property 1 or Property 2. 

For each i, hi(w) > for each w in C such that Wi(0) = 0, u>i(s) > for 
s < 0, and Wj(s) - for s < andj ¥> i. 
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